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ABSTRACT 

We demonstrate the Sunyaev-Zel'dovich (SZ) effect imaging capabilities of the Combined Array for 
Research in Millimeter- wave Astronomy (CARMA) by presenting an SZ map of the galaxy cluster 
RX J1347.5— 1145. By combining data from multiple CARMA bands and configurations, we are able 
to capture the structure of this cluster over a wide range of angular scales, from its bulk properties to 
its core morphology. We find that roughly 9% of this cluster's thermal energy is associated with sub- 
arcminute-scale structure imparted by a merger, illustrating the value of high-resolution SZ measure- 
ments for pursuing cluster astrophysics and for understanding the scatter in SZ scaling relations. We 
also find that the cluster's SZ signal is lower in amplitude than suggested by a spherically-symmetric 
model derived from X-ray data, consistent with compression along the line of sight relative to the 
plane of the sky. Finally, we discuss the impact of upgrades currently in progress that will further 
enhance CARMA's power as an SZ imaging instrument. 
Subject headings: galaxies: clusters: general 



1. INTRODUCTION 

Galaxy clusters are the largest gravitation ally bound 
systems in the universe, and have taken nearly a Hubble 
time to form. They therefore have the potential to act 
as powerful probes of cosmology if systematic errors can 
be controlled. Precision cluster cosmology will require a 
deep understanding of cluster astrophysics, particularly 
as it relates to the hot gas of the intracluster medium 
(ICM). The most detailed studies of the ICM have thus 
far been performed by X-ray telescopes, which are sensi- 
tive to the bremsstrahlung emission from the 10 7 -10 8 K 
gas. The Sunyaev-Zel'dovich (SZ) effect is a complemen- 
tary probe of the ICM. The amplitude of the SZ signal 
depends on the line-of-sight integral of n e T, while the 
X-ray surface brightness depends on n^, so sensitive SZ 
measurements can access tenuous gas outside the cluster 
core and directly measure pressure disturbances. Fea- 
tures found commonly in the outer regions of clusters, 
such as shocked gas from mergers, may therefore be eas- 
ier to detect using the SZ effect than using X-rays. More- 
over, the combination of X-ray and SZ data can be used 
to obtain a more complete picture of the ICM thermo- 
dynamics. 

To take advantage of these opportunities, advances in 
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SZ imaging capabilities are needed. Measurements of 
the SZ effect have become routine over the last decade, 
but the full potential of the SZ effect as a probe of clus- 
ter physics remains largely unexploited due to technical 
challenges: the combination of high sensitivity and large 
angular dynamic range required for detailed SZ imaging 
has proven difficult to achieve with existing instruments. 
As a result, the use of the SZ effect has been limited pri- 
marily to studies where resolved imaging is unnecessary. 

The small number of higher-resolution SZ images ob- 
tained to date have served to demonstrate the utility of 
the technique . However, sin g le- dish measurements such 
as those by Korngut et al.| ( |2Q11| ) can suffer from ra- 
dio point source contamination and have been limited 
to scales < 45" by the necessity of filtering out modes 
contaminated by atmospheric noise. Multi- dish SZ mea- 
surem ents by arrays such as ATCA (e.g., Malu et al. 
2010) can constrain and remove point sources using the 
inherent spatial filtering ability of interferometers, but 
most millimeter-wave interferometers lack sensitivity at 
arcminute angular scales where the SZ cluster signal is 
largest. 

CARMA is a heterogeneous interferometric array con- 
sisting of 23 antennas with diameters of 3.5, 6.1, and 
10.4 m operating at 1 cm, 3 mm, and 1 mm. This partic- 
ular combination of antennas and bands makes CARMA 
a uniquely powerful SZ instrument: its 3.5 m antennas 
can be placed in a compact configuration sensitive to 
arcminute-scale emission, and its 6.1 and 10.4 m dishes 
can be used to obtain the sensitivity necessary to resolve 
smaller angular scale SZ features. In this work, we make 
use of CARMA data from three array configurations and 
two bands to obtain an SZ image of the galaxy cluster 
RX J1347.5-1145. These data represent the highest- 
fidelity picture of a galaxy cluster ever obtained using 
the SZ effect. 

This paper is organized as follows: Section [2] provides 
background on RX J1347.5— 1145, Section[3]describes the 
observations and data reduction, and Section [4] discusses 
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the modeling and deconvolution method. We present 
our results and compare them to previous measurements 
in Section [5| and review the conclusions and discuss 
prospects for future work in Section [6j 

2. RX J1347.5-1145 

The target for these observations is the cluster 
RX J1347.5— 1145, an object that has been character- 
ized extensively using a variety of tec hniques. First dis - 
covered by the ROSAT all-sky survey flVoges et al|l 999 ) , 
RX J 1347. 5— 1145 is the most luminous cluster known in 
the X-ray sky, and has been meas ured by several X -ray 
instruments inc luding Chandra (| Allen et al. 2002) and 
XMM-Newton flGitti ~fe Schindler||2004[ ). Optical obser- 
vations have revealed the presence of two cD galaxies, 
one coincident with the X-ray emission peak and one di- 
rectly to the east. The system has also been found to host 
a radio mini-halo (Gitti et al. 2007). Its gravitational 
potential has been probed using both strong and weak 
gravitationa l lensing (e.g., |Miranda et al.|[2008| |Bradac| 
|et al.||2008| ). These mult i- wavelength observations indi- 
cate thaTRX J1347.5-1145 is a massive (> 10 15 M ) 
cluster at redshift z = 0.4510 which has recently under- 
gone a merger with a smaller object. The cluster's SZ sig- 
nal has been measured using single-dish (Kom atsu et al. 
2000||Kitayama et al.||2004l|Korngut et aH[5Dll|) and in- 



terlerometr ic (Carlstrom et al. ' 2002[ |Bonamente et 
20081 12011|) imaging instruments, and its spectrum n 



aL 
spectr um near 



the thermal SZ null has also been characterized (Zemcov 
et al. 2012| ). The higher angular resolution measurements 
have revealed a compact region of very hot (~ 20 keV) 
gas to the southeast of the X-ray emission peak, while the 
low-resolution data indicate a smaller arcminute-scale SZ 
signal than suggested by a spherical fit to the X-ray data. 
The existence of a central radio-bright AGN, along 
with the limited angular dynamic range of most SZ in- 
struments, has complicated efforts to bring SZ data to 
bear on understanding this system. The CARMA data 
we present help to overcome both limitations. 

3. OBSERVATIONS AND DATA REDUCTION 

RX J1347.5— 1145 was observed with three different 
sets of CARMA antennas at two wavelengths: an 8- 
element sub-array consisting of 3.5 m antennas at 1 cm 
("CARMA-8"), a 15-element sub-array consisting of 
6.1 m and 10.4 m antennas at 3 mm ("CARMA- 15"), 
and the full 23-element array at 3 mm ("CARMA-23"). 

The CARMA-8 data were obtained in August 2009 
and totaled 25.7 hours of unflagged, on-source time. 
The center frequency was 31 GHz with a bandwidth of 
8 GHz, and the target R.A. and decl. were 13:47:30.7 and 
-11:45:08.6 in J2000 coordinates. The array was config- 
ured with six elements in a compact array sensitive to 
arcminute-scale SZ signals, and two outlying elements 
providing simultaneous discrimination for compact radio 
source emission. The compact array and longer base- 
lines sample w-spacings of 350—1300 A and 2—7.5 kA, 
respectively. The data were reduced using the Sunyaev- 
Zerdovich A rray (SZA) pipeline described in Muchovej 
eTaLl([2007|. 



The CARMA- 15 data were obtained in July 2010 and 
totaled 7.4 hours of unflagged, on-source time. The 
center frequency was 90 GHz with a bandwidth of 
8 GHz, and the target R.A. and decl. were 13:47:32.0 



and -11:45:42.0 in J2000 coordinates. The CARMA-15 
array was pointed slightly southeast of the CARMA-8 
phase center, directly toward the region of hot gas. The 
antennas were in the E configuration, the most compact 
standard positions for the 6.1 and 10.4 m ant ennas. Data 
reduc tion was performed using MIRIAD (|Sault et al.| 



1995]). 

We obtained C ARMA-23 data as part of a commission- 
ing run in February 2011. All 23 antennas were operated 
at a center frequency of 86 GHz and were attached to 
the CARMA spectral line correlator, which operated at 
a reduced bandwidth — 4 GHz for the double-sideband 
receivers on the 10.4 and 6.1 m antennas, and 2 GHz 
for the single-sideband receivers on the 3.5 m antennas. 
A total of 8.6 hours of unflagged on-source data were 
obtained; a relatively high fraction of the 3.5 m data 
were flagged due to hardware issues in the commission- 
ing run that have since been corrected. The target R.A. 
and decl. were the same as for the CARMA-8 data, and 
the array configuration was approximately the combined 
CARMA-8 and CARMA-15 configurations. As with the 
CARMA-15 data, we reduced the C ARMA-23 data using 
MIRIAD. 

For the CARMA-15 and CARMA-23 arrays, we treat 
each baseline type (10. 4mx 10.4m, 10. 4mx6.ini, etc.) 
separately to properly account for the differing pri- 
mary beams. We therefore have ten data sets: one for 
CARMA-8, three for CAMRA-15, and six for CARMA- 
23. We apply a cutoff in uv radius for each data set, 
using the data beyond the cutoff only to constrain the 
point source emission. The cutoff is chosen to exclude 
portions of the uv plane which are poorly sampled for a 
given baseline type. 

The outputs of the data reduction pipelines consist 
of flagged, calibrated visibilities V( u,v). All ab solute 
flux calibration is performed using the |Rudy|Ql987| ) Mars 
model, and is accurate to 5%. To combine data from dif- 
ferent bands, we define y(u,v), the Fo urier- space coun- 
terpa rt to the Compton y parameter (Carlstrom et al. 
[2002] ): 

V v (u,v) 



y{u,v) 



9{v,(Te))Io 



(1) 



where (T e ) = 12 keV is the mean ICM electron tem- 
perature of the cluster 9 . Our ten y(u,v) data sets are 
summarized in Table [I] 

The data weight distribution in the uv plane for the 
combined data sets is shown in Figure [I] To illustrate 
the contributions of each of the three sub-arrays, we show 
the measured y(u, v) binned in uv radius in Figure [2] 
CARMA-8 at 1 cm constrains the cluster at small uv 
radius (large angular scale) where the signal is largest, 
while CARMA-15 and CARMA-23 at 3 mm are sensitive 
to the large uv radius (small angular scale) substructure 
of the cluster. 

4. MODELING AND DECONVOLUTION 

We wish to combine the information in all ten sets of 
visibility data to form a single image of the cluster. The 
first step in this process is to remove the radio point 

9 We make this approximation due to our incomplete of knowl- 
edge of the true T e (x, y, z). The relative value of g[y, T e ) between 
1 cm and 3 mm is nearly independent of T e , varying by just 2% 
between 5 and 15 keV. 
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Fig. 1. — Normalized data weight distribution in the uv plane for the union of data sets described in Table pH Weights are calculated 
from the inverse variance of each visibility, scaled by the SZ intensity spectrum. The iw-plane extent of each visibility weight is determined 
from the cross-correlation of the illumination patterns of the corresponding antennas, providing a more complete view of the uv sampling 
in the heterogeneous array. The weights are well-matched to the cluster signal — which is largest at small uv radius — except for a relatively 
under-sampled region around ~ 2 kA. This region of the uv plane will be well-measured by the 23-element CARMA array at 1 cm, which 
is currently under development. 
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Fig. 2.— Radially binned visibilities from CARMA-8 at 1 cm, CARMA-15 at 3 mm, and CARMA-23 at 3 mm. The CARMA-8 data 
measure the cluster signal at large angular scales, while the other sub-arrays measure the smaller substructure of the cluster. 
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source emission, which is accomplished by fitting a point 
source model to the visibilities above the uv cutoff. Since 
radio source fluxes often vary considerably with time, 
identical fluxes should not be expected in the CARMA-15 
and CARMA-23 despite their similar central frequencies. 
We therefore demand consistency within each individ- 
ual sub-array, but allow the source flux to vary between 
sub-arrays. We find a best-fit centroid consistent with a 



within the central region with an a = —2.331 power law 



source in the NVSS catalog (Condon et al. 1998) at an 



R.A. and decl. of 13:47:30.7 and -11:45:08.6. The flux 
of this source was found to be 8.9±0.5 mJy in the 1 cm 
CARMA-8 data, 4.2±0.2 mJy in the 3 mm CARMA-23 
data, and 4.0±0.2 mJy in the 3 mm CARMA-15 data. 
We subtract the best-fit point source models from the 
visibility data sets before proceeding. 

We next apply a phase shift to the CARMA-15 visi- 
bilities to establish a single phase center for all ten data 
sets. At this stage, all of the visibilities can be combined 
to form a single dirty map (Figure [3]). However, this map 
is difficult to interpret because it includes data with dra- 
matically different primary beams, and because noise at 
small angular scales masks the arcminute-scale emission. 
In order to proceed, we must use some image deconvo- 
lution method such as the Maximum Entropy Method 
(MEM) or CLEAN. For simplicity, we choose to build a 
model of the cluster iteratively using the CLEAN algo- 
rithm, though we note that MEM also holds promise for 
future work. 

Our algorithm successively builds up a single model 
of the cluster in image space by CLEANing individual 
data sets. We begin by making a dirty map from the 
first data s et y\{u, v). We then run the Hogbom CLEAN 



i 



algorithm (Hogbom 1974b on this map with gain 0.05, 



stopping at 2.5cr. The CLEAN components are corrected 
for the appropriate primary beam, and a restored map is 
generated. This restored map (in units of Compton y per 
pixel) is our initial cluster model. The cluster model is 
then multiplied by the primary beam of the next data set, 
y2(u,v), and the result is Fourier transformed and sub- 
tracted from 2/2 (^? v). From these model-subtracted visi- 
bilities, a dirty map is produced and Hogbom CLEANed. 
The CLEAN components are corrected for the primary 
beam and restored, and the restored map is added to the 
cluster model. We repeat this process on all ten data 
sets, resulting in a cluster model that incorporates infor- 
mation from all of our data. We continue refining the 
model by making additional passes through all ten data 
sets until each model-subtracted dirty map is consistent 
with noise. Three iterations are found to be sufficient. 
Our final cluster model is then added to the weighted 
average of the ten model-subtracted dirty maps to form 
a CLEAN map (Figure EL Though the resolution of the 
map is not well-defined, the smallest restoring beam of 
the ten used to construct the model is 10'/ 6 x 16'/9. 

To assess the accuracy and flux recovery of our de- 
convolution technique, we construct a heuristic model of 
RX J 1347. 5— 1145. We generate mock data using this 
model, repeat the procedure described above on the sim- 
ulated visibilities, and compare the results to the input 
model. Modeling the SZ signal requires knowledge of 
the density n e and the temperature T e of the ICM, both 



of which ca n be approximated using the results of | Allen 



et al. ( 2QQ2J) (hereafter [A02[ ) derived from Chan dra X-ray 
data. For the electron density, we use the |A02| density fit 



in the outer region inferred from the A02 surface bright- 
ness profile fit. We piece together the inner and outer 
regions by requiring that n e (r) be co ntinuous. For the 
electron temperature, we use the A02 measured temper- 
ature profile out to the outer radius of their largest bin, 
and use the average temperature of 12.0 keV at larger ra- 
dius. We assign a temperature of 18.0 keV in the shocked 
region of the southeast quadrant. Finally, we truncate 
both the temperature and density profiles at r2oo- 

Using these approximations for the electron density 
and temperature, we integrate along the line of sight 
to produce a simulated Compton y map, and Fourier 
transform to produce simulated visibilities y S im(u, v). We 
then randomize the phases in our visibility data y(u, v) 
to produce simulated noise y n oise( u , v ) an d run our iter- 
ative CLEAN algorithm on y S im(u, v) + y n oise(u, v). We 
find that our algorithm accurately reproduces the mor- 
phology of the input model, and that it recovers 69%, 
78%, and 92% of the integrated flux within r2oo? r 500? 
and r25oo, where th e ta values are determined from the 
Allen et al.| Q2002D NFW model fit. The flux recovery 
ratio at larger radii is highly sensitive to the assumed 
electron density power law index. This is due to the fact 
that r2oo fc> r this cluster corresponds to an angular scale 
of 5.75, comparable to the 10.7 FWHM 3.5 m primary 
beam at 1 cm. 

5. RESULTS 

The map shown in Figure [4] can be understood as a 
relaxed cluster SZ signal with additional sub-arcminute 
structure in the southeast quadrant imparted by the 
merger event. The imaging of both the extended and 
compact structure in our SZ map is of significantly higher 
fidelity than previous measurements due to the ability 
of CARMA to remove the central point source and to 
the large angular dynamic range of the combined arrays. 



In contrast with pr evious work (Mason et al.|[2010 



Ko- 



matsu et al. 2000), we find that the peaks of the SZ 



and X-ray signals are coincident. Figure [6] shows a com- 
parison between the CARMA SZ map and the Chandra 
X-ray pressure map ( Bradac et al.||2008J ). Though a full 
multi-wavelength reconstruction is beyond the scope of 
this paper, it is obvious that the two techniques produce 
maps with consistent morphologies. 

To separate the relaxed and disturbed components, we 
make use of the p ressure profi le fit to the Chandra X-ray 



data described in Allen et al. 



(2008), in which the south- 



east quadrant was excised. We first project the pressure 
profile along the line of sight to produce an integrated 
pressure map. For each sub- array, we multiply this map 
by the appropriate primary beam, convert to Compton 
?/, Fourier transform, and subtract the result from the 
visibility data. We then iteratively build a model of the 
remaining SZ signal, following the procedure described 
in Section [4j We allow the X-ray pressure profile to 
be scaled by a multiplicative constant to compensate 
for cluster projection effects and calibration errors. The 
scale factor is chosen so as to produce no net CARMA- 
8 (arcminute-scale) signal in the iteratively-determined 
model. 

Removing this estimate of the relaxed signal from our 
visibility data allows us to focus on the sub-arcminute- 
scale signal resulting from the merger event. The result 
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Fig. 3. Dirty map of RX J1347.5-1145 using "robust" 
The map is in units of Compton y x 10 13 /beam. 

is shown in Figure [5] The total Compton y signal re- 
covered in this map, which corresponds to the fraction 
of the thermal energy in the cluster ICM associated with 
the merger-related substructure, is 9.1% of the total re- 
covered from the map in Figure [4j 

The value of the multiplicitive constant by which the 
X-ray pressure is scaled can also yield information about 
the characteristics of the cluster. Since the scale factor is 
determined by requiring consistency with the CARMA-8 
data, its value depends upon the three-dimensional mor- 
phology of the ICM on large angular scales. Due to the 
different dependencies of the X-ray surface brightness 
and the SZ signal on density, a cluster more (less) ex- 
tended along the line of sight than in the plane of the 
sky would yiel d a scale factor of g reater than (less than) 



visibility weighting. The synthesized beam is shown in grey in the lower left. 



one (see e.g., Grego et al 



find a scale factor of 0.58 



2004). For this cluster, we 
indicating a cluster which is 
strongly compressed alon g the line-of-sight directi on. A 



similar ratio is reported in |Bonamente et aL (|2011 ) using 
the C ARMA-8 data reported here. [Chakrabarty et al.| 
(|2008| ), using a combination of X-ray and SZ data, also 
find that RX J 1347.5— 1145 is compressed along the line 
of sight with an axis ratio of ~ 5. A simple comparison 
of the X-ray surface brightness profile and our data, as- 
suming an isothermal cluster, implies compression along 
the line of sight by a factor of roughly three-to-one — 



a fairly extreme value, t hough less than suggested by 
Chakrabarty et al.| ( |2008 ). 

Clumping of the ICM, i.e., a systematic discrepancy 
between (n^(r)) and (n e (r)) 2 , could also lead to differ- 
ences between the SZ and X-ray signals. Clumping is 



observed in the outskirts of simulated clusters (Nagai & 



Lau 2011), and is implied by observations of flattened 



entropy profiles and gas mass fractions apparently in ex- 
cess of the cosmic mean in some clusters (e.g.,|Simionescu 



!01ip . 
eclus 



let al.|2011|). However, in both cases the clumping occurs 



at large cluster radii (at least > rsoo), and is thus unlikely 
to affect the normalization of the X-ray model used here 
(which was fit to data at r < r^oo) at the level required 
to explain the offset. 

The difference between the CARMA and Chandra 
pressure estimates could also arise in principle from cali- 
bration or systematic errors in the SZ or X-ray data. As 
a cross-check of our CARMA calibration, we compared 
o ur data to previou s BIMA SZ measurements reported 
in [Bonament e et aL| ( 2008 ) , finding consistent integrated 
Y values and binned visibility data. Since the Chandra 
calibration is unlikely to be mistaken at this level, and 
since major mergers such as RX J1347.5— 1145 are not 
expected to be spherically symmetric, we suggest that 
line-of-sight compression is most likely to be the domi- 
nant effect. However, ICM clumping and calibration er- 
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Fig. 4. CLEAN map of RX J1347.5-1145. The map is in units of Compton y x 10 15 /pixel, where the pixel size is 0.5" x 0.5". The 
smallest CLEAN beam used to construct the model is 10". 6 x 16". 9 
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Fig. 5. — CLEAN map of RX J1347.5— 1145 with the scaled relaxed X-ray pressure profile subtracted from the visibilities. The map is 
in units of Compton y x 10 15 /pixel, where the pixel size is 0.5" x 0.5". 
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rors may also be contributing to the discrepancy; a more 
complete explanation will require a joint analysis of the 
two data sets. 

6. CONCLUSIONS AND FUTURE WORK 

We have demonstrated the ability of CARMA to mea- 
sure the SZ signature of galaxy clusters at high sensitivity 
across a wide range of angular scales. By combining data 
from three CARMA configurations and two frequency 
bands, we measure the arcminute-scale SZ signal as well 
as the pressure substructure of RX J1347.5— 1145. The 
large angular dynamic range of the CARMA data and the 
capability to constrain and remove point sources make 
our measurement a significant improvement in ima ge fi- 
delity over previous work (e.g., |Korngut et al.|2011[). By 
comparing o ur data to the X-ray measurements of |Allen| 



et al. 



(120081), we are able to determine that ~ 9% of the 
SZ signal is localized in the disturbed region of the clus- 
ter, and that the system is likely compressed along the 
line of sight relative to the plane of the sky. 

Although our results demonstrate that CARMA in its 
current form is a highly capable SZ instrument, several 
upgrades will soon bring about significant enhancements. 
An effort is currently underway to equip all antennas 
with 1 cm receivers and expand the correlator bandwidth 
to 8 GHz for a 23-element array. The 3 mm CARMA- 
23 array described in this work provided sensitivity from 
uv radii of ~ 2.0 to lOkA; the upgraded array placed 
in the same configuration but operated at 1 cm will pro- 
vide higher- sensitivity coverage from ~ 0.35 to 3. 3k A. As 
shown in Figure |2j this corresponds to the portion of the 
uv plane where the SZ signal is large. A planned upgrade 
to more sensitive 3 mm receivers will allow the SZ signal 
at finer angular scales to be measured more precisely. 

The sensitivity to smaller angular scale SZ structures 
provided CARMA's larger telescopes can be directed to- 



ward regions of interest, as was done in this study by 
pointing them toward the region of hot gas to the south- 
east of the cluster center. With the increased sensitivity 
enabled by the ongoing upgrades, this technique can be 
used to search for shock-enhanced features in the out- 
skirts of clusters. Mosaicking can also be used to pro- 
vide sensitivity to small-scale structures over the entire 
cluster. 

Taken together, these upgrades and observing strate- 
gies will allow CARMA to image clusters precisely and 
efficiently over a wide angular dynamic range, making it 
possible to fully exploit the power of the SZ effect as a 
probe of cluster astrophysics and precision cosmology. 
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TABLE 1 
CARMA Observations of RX J1347.5- 



1145. 



Array 


Ant 1 


Ant 2 


Frequency 


uv cutoff 


Noise 


Minor axis 


Major axis 


Beam P.A. 








(GHz) 


(kA) 


mJy/beam 


arcsec 


arcsec 


Degrees 


CARMA-8 


3.5m 


3.5m 


31 


2000 


0.16 


100.1 


128.8 


-33.4 


CARMA-15 


10.4m 


10.4m 


90 


10000 


0.38 


13.4 


16.8 


72.9 


CARMA-15 


10.4m 


6.1m 


90 


10000 


0.19 


12.8 


16.5 


-44.8 


CARMA-15 


6.1m 


6.1m 


90 


10000 


0.30 


15.1 


37.4 
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86 


10000 
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88 
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0.80 
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18.3 


59.1 
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88 


10000 


0.76 


10.6 


16.9 


-81.1 
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2.64 


31.4 


67.5 


-44.3 
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Fig. 6. — The C LEANed CAR MA SZ map (left) cropped and scaled for direct comparison with the published X-ray-derived integrated 
pressure map from Brad ac et al.||2008| (right). 



